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Abstract 



The temperature dependence of the moment of inertia of a dilute Bose gas 
confined in a harmonic trap is determined. Deviations from the rigid value, 
due to the occurrence of Bose-Einstein condensation, reveal the superfluid 
behaviour of the system. In the noninteracting gas these deviations become 
important at temperatures of the order of TcN^^/^"^ . The role of interactions 
is also discussed. 
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Evidence of Bose-Einstein condensation in magnetically trapped gases of ^"^Rb and 
^Li atoms has been recently reported. These measurements are expected to open further 
stimulating perspectives in the study of the phenomenon of Bose-Einstein condensation 
(BEC) both from the experimental and theoretical point of view [^. In this context the 
understanding of the superfluid behaviour of a Bose condensed atom cloud is a question of 
primary importance that we plan to investigate in the present letter. 

A natural way to discuss superfluidity in a confined system is to focus on its rotational 
properties [Q. For a macroscopic system the moment of inertia is given by the rigid value 
unless it exhibits superfluidity. Crucial deviations from the rigid motion occur in rotating 
liquid helium below the lambda temperature [|5[. Also finite systems, like some deformed 
atomic nuclei p and helium clusters are known to exhibit important superfluid effects 
in the moment of inertia. 

The moment of inertia 0, relative to the z-axis, can be defined as the linear response of 
the system to a rotational field Hext = —^Jz, according to the formula 

< >= cue (1) 

where Jz = J2i{xiPi —yipf) is the ^-component of the the angular momentum and the average 
is taken on the state perturbed by Hext- For a classical system the moment of inertia takes 
the rigid value 

Grig = mN <x^ + > (2) 

where is the number of atoms in the trap. Vice versa the quantum mechanical determi- 
nation of O is less trivial. It involves a calculation of dynamical properties of the system 
and, according to perturbation theory, can be written as 

= ^ \^ „-Em/kBT \< ml Jz\n>\'^ , . 

Z E - E 

^ n,m 

where m (n) and Em [En) are eigenstates and eigenvalues of the unperturbed hamiltonian 
and Z is the partition function. Generalization of (|^) to the grand canonical ensemble is 
straightforward. 
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In the following we calculate the moment of inertia of a dilute Bose gas confined in 
a harmonic trap. We first consider the simplest case of an ideal gas described by the 
hamiltonian 

^ = E ^ + E + -ii + ■ (4) 

i i 

The role of interparticle interactions will be discussed in the second part of the work. For the 
model hamiltonian @ the moment of inertia is explicitly evaluated by solving the equation 

[H,X] = J, (5) 

for the operator X which, according to (0), determines the moment of inertia through the 
relation 6 =< >. The explicit form of the operator X is found to be 

X = - J + u^Dx^y. + ^p!] (6) 

and the moment of inertia takes the form 

= ^^[(< y'>-<x' >Kcol + col) + 2^ <y"> < >)] (7) 

where the average is taken at statistical equilibrium in the absence of the perturbation Hf.xt 
Notice that expression makes sense only if uo^ 7^ cOy. 

Equation (|^) is an exact result for the hamiltonian @). It applies to any value of and 
holds both for Bose and Fermi [0] statistics. For the explicit evaluation of the moment of 
inertia in a Bose gas it is useful to evaluate the average radii < > and < > in the 
so called macroscopic limit where one separates the contribution arising from the conden- 
sate from the one of the other excited states, whose quantum numbers are approximated 
by continuous variables (semiclassical approximation). This description is appropriate for 
temperatures fc^T much larger than the oscillator energies huj^i huy, and yields the 



following result for the depletion of the condensate |11 



N - NoiT) = (8) 



with the critical temperature given by 



Here Uq = u^UyUz and Q{ri) = s"^/ (exp s — l)ds. For the square radius < > one finds 
the result 

and analogously for < >. The first term of this equation gives the contribution to < > 
arising from the particles in the condensate and scales as l/uj^- The second one is the 
contribution from the non condensed particles and scales as I/ujI- When inserted into 
the contributions to < a;^ > and < > which scale as l/uj^ and l/cUy, give rise to the 
irrotational value for the moment of inertia. Conversely contributions which scale as l/uj"^ 
and l/ujy yield the rigid value. The final result is 

e = elm<x^ + >o iVo(T) + m < + (AT - iVo(T)) (11) 

where the indices <>o and <>nc mean average on the Bose condensed and non condensed 
components of the system. The quantity 

.„ = (12) 

is the deformation parameter of the Bose condensed density. In terms of the oscillator 
frequencies, it is given by eo = {yjy — ijJx)l{y>y + '^x)- 

The physics contained in equation (|ll|) is very clear. The moment of inertia is the 
sum of two terms arising from the atoms in the condensate (superfluid component), which 
contribute with their irrotational flow, and from the particles out of the condensate (normal 
component) which rotate in a rigid way. These two distinct contributions are at the origin 
of an interesting T-dependence of the moment of inertia. In fact above T^, where Nq{T) = 0, 
the moment of inertia takes the classical rigid value (^, while at T = 0, where all the atoms 
are in the Bose condensate, is given by the irrotational value Qirrot = e^Qrig- The explicit 
T dependence is easily calculated using the harmonic oscillator result 



<x^ + y^ >nc ^ 2 g(3) u;l + LJ^y ^^^^ 
< x2 + ?/2 >p 3 Q(2) + ujy)uJxUJy ^ 

Assuming c^a. ~ 7^ uj^ and taking the values Q{2) = 2.4 and (5(3) = 6.5, one obtains 
the expression 



e„g 1 - (T/T,)3 + 1.7(T/T,)4(AAr)i/3 

where A = uOz/oJ^- Near the moment of inertia deviates very little from the rigid value not 
only because of the smallness of the number of atoms in the condensate, but also because of 
the large A^^^^ factor in the denominator of Eg . (|14|) . This is the consequence of the fact that 
the radius of the atoms in the condensate is much smaller than the one of the atoms out 
of the condensate. In order to observe superfluid effects in the moment of inertia one has 
to go to smaller temperatures of the order of T ~ T^N^^^^"^ where the term proportional to 
N^/"^ in the denominator of Eq.([T^) is of the order of unity. It is worth noting that, even for 
large values of A^, the relevant regime of temperatures is not much smaller than Tc and can 
still be large compared to the oscillator temperature. In the figure we report the prediction 
of Eq.([T^) for two different values of A^ and A = (this value characterizes the magnetic 
trap of 0). 

Result (^) for the moment of inertia has been derived for an ideal Bose gas confined 
in a harmonic trap. The same result is expected to hold also for a weakly interacting gas 



where the mean field description permits to obtain coupled equations for the wave 

functions of the condensate and of the single particle excited states. The irrotational flow for 
the condensate then follows from general properties of the equation for the phase. On the 
other hand the rigid flow for the particles out of the condensate follows from the use of the 
semiclassical approximation which should be applicable in a useful range of temperatures 
1^ . Interactions can neverthelss modify the value of the moment of inertia (|Tl|) by changing 
the temperature dependence of Nq(T) as well as the value of the square radii < + y"^ > 



mm- 

In order to give a first estimate of these effects we have considered the Hartree-Fock 
equations at temperatures such that the thermal depletion of the condensate is small and 



one can consequently ignore the contribution of excited atoms to the self-consistent effective 
field. The temperature should be however large enough in order to apply the semiclassical 
approximation to the excited states. This regime is the most interesting for the investigation 
of superfluid effects in the moment of inertia. Furthermore, in order to obtain an explicit 
estimate, we have considered the case of strongly repulsive interactions where the the kinetic 



term in the Gross-Pitaevskii equations for the condensate can be neglected ||12|,|T3 . 

Under the above conditions (small thermal depletion, semiclassical approximation for 
the excited states and strongly repulsive interactions) the Hartree-Fock equations at finite 
temperature yield the following simplified expressions for the condensate (po) and non con- 
densate (pnc) densities: 

in 

po{r) = 2 - VextWifi - Vext) (15) 

47r/i a 

and 

Pnc(r) = ( ^^^f )^^^^3/2(exp(- I p - Vext \ /ksT)) . (16) 

In the above equations a is the scattering length, 6 is the unit-step function, g^i2 is the usual 
Bose function, Vext = ^iy^lx^ + ^yV"^ + ^1^^) is the harmonic potential and 

/i = hlh^uolNh^afl^ (17) 



is the T = value of the chemical potential, fixed by the normalization of po- From Eqs.(15) 
and ( |T^ ) one gets the useful results (valid if (A^ — A^o) « N): 

N-No _T' p 

- Yf'^k^^ ^^^^ 

and 

<x^ + y^>o 3Fo(^) p ^ ^ 

with Fnifi/ksT) = J,°° dss^+^gs/2{exp{- \ - fi/kBT \)/ J,^ dss^g3/2{exp{-s^)). 

In Eq.(|TB|) Tc is the critical temperature of the noninteracting Bose gas. We assume 
that this value is not significantly affected by the interaction. 
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With respect to the predictions (|) and (|T^) of the ideal Bose gas, the interactions 
increase the thermal depletion of the condensate {N — Nq{T)) due to the (positive) chemical 
potential entering Eq.(|l^). On the other hand they can significantly decrease the ratio 
< + >nc I < + y"^ >o since the radius of the condensate density is very sensitive to 
the repulsive forces. As an example we have chosen N = 20000 and the values a = IOOoq {clq 



is the Bohr radius), ax = y = 1-65 x 10~^cm ~ ay and A = y/S, typical of the rubidium 
trapped gas of Ref. With these values we find ksTc = 36.1huJx and ^ = ll.GhuJx. At 
T = Tc/3 the thermal depletion is increased from 0.04 (prediction of the noninteracting 
model) to 0.09, the ratio < + y'^ >nc / < x"^ + y"^ >o is decreased from 21.7 to 4.0 and the 
ratio Q/Qrig from 0.45 to 0.3, thereby showing that in this case the effects of superfluidity 
in the moment of inertia are reinforced by the interactions. 

A more systematic investigation in a wider range of temperatures and values of the scat- 
tering length (including the interesting case of attractive forces) requires the self-consistent 
solution of the Hartree-Fock equations at finite temperature. The present analysis, however, 
already permits to draw a first important conclusion, i.e. that superfiuid effects in the mo- 
ment of inertia of a magnetically trapped Bose gas should be observable at temperatures 
not dramatically smaller than the critical temperature for Bose-Einstein condensation. 

A final question concerns the onset of rotational instabilities occurring at large frequencies 
uj. For an axially symmetric gas the stability is ensured below the critical frequency = 
min{ej^/ Jz) where ej^ is the energy of a general excitation carrying angular momentum J^. 
In the noninteracting case the value of ujc coincides with the oscillator frequency uj^- In the 
presence of interactions one should calculate the dispersion law of the elementary excitations 
as well as the energy needed to create a vortex line [0. Work in this direction is in progress. 

Useful discussions with D. Brink, E. Cornell and F. Laloe are acknowledged. 
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FIGURES 

FIG. 1. Temperature dependence of the condensate fraction (full line) and of the moment of 

inertia for N = 2000 (dashed) and N = 20000 (dot-dashed) atoms in the noninteracting gas (see 
the text). 
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